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SUMMARY 

We report palaeointensity estimates obtained from three Ohgocene volcanic sections 
from the Kerguelen Archipelago (Mont des Ruches, Mont des Tempetes, and Mont 
Rabouillere) . Of 402 available samples, 102 were suitable for a palaeofield strength de- 
termination after a preliminary selection, among which 49 provide a reliable estimate. 
Application of strict a posteriori criteria make us confident about the quality of the 
12 new mean-flow determinations, which are the first reliable data available for the 
Kerguelen Archipelago. The Virtual Dipole Moments (VDM) calculated for these flows 
vary from 2.78 to 9.47 10^^ Am^ with an arithmetic mean value of 6. 15 ±2.1 10^^ Am^. 
Compilation of these results with a selection of the 2002 updated lAGA palaeointensity 
database lead to a higher (5.4±2.3 l^)"^"^ Arn^) Oligocene mean VDM than previously re- 
ported dGoguitchaichvih et al., 20011 [Riisager, 1999D , identical to the 5.5±2.4 lO^^Am^ 



mean VDM obtained for the 0.3-5 Ma time window. However, these Kerguelen palaeoin- 
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tensity estimates represent half of the rehable Ohgocene determinations and thus a bias 
toward higher values. Nonetheless, the new estimates reported here strengthen the con- 
clusion that the recent geomagnetic field strength is anomalously high compared to that 
older than 0.3 Ma. 
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1 INTRODUCTION 

Numerous studies have been carried out to increase our knowledge of geodynamo physics, 
but even so we still do not know the detailed mechanism of the generation of the Earth's 
magnetic field, and even less about the processes that produce secular variation, excursions 
and reversals (e.g., Jacobs ()1994|) : Merrill & McFadden ()1999|l ). To better understand the 
geomagnetic field we need to be able to go back in time in order to observe its changes and 
to obtain long-term global characteristics. This is possible with some rocks which recorded 
the Earth's palaeomagnetic field during their formation. Volcanic rocks, in particular, fur- 
nish a global knowledge of the geomagnetic field because they contain information on both 
the direction (inclination and declination) and the strength of the palaeofield. However, the 
most reliable methods for absolute palaeointensity determination, Thellier & Thellier ()1959jl 
and its modified version proposed by Coe ()1967a|) . are time consuming because of the strict 
conditions which have to be checked to validate the determinations. Moreover, many vol- 
canic rocks turn out to be unsuitable for palaeointensity determination. For these reasons, 
reliable palaeointensity data are difficult to obtain and are particularly rare. Only 1.5 de- 
terminations per million years between and 300 Ma ( Selkin fc Tauxe, 2000D are available 



when combining the updated lAGA 1999 data-set and the Scripps submarine basaltic glass 
databases. It is obvious that more palaeointensity data are needed for a better understand- 
ing of the geomagnetic dynamo. 

* Corresponding author. Tel: +33 467 14 39 38; Fax: +33 467 14 36 03. e-mail address: camps@dstu.univ-montp2.fr 
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This study on three basaltic sections sampled in the Kerguelen Archipelago (49.9°S, 70°E) 
aims to estimate more accurately the palaeomagnetic strength of the geomagnetic field in 
the 28-30 Ma time interval recorded by these lavas. It thus complements the palaeomagnetic 
directions recently published on the same sections ( Plenier et al., 2002D . By combining these 



new determinations with selected records issued from preexisting palaeomagnetic databases 
(Tanaka et al. (fT?M|l . updated by Perrin & Shcherbakov (nW7|l and Perrin et al. (ITO^ l 
we will propose a more robust estimation of the Oligocene palaeofield strength and we will 
discuss the long-term characteristics of the geodynamo. 

2 GEOLOGY AND SAMPLING 

The Kerguelen Archipelago lies on the northern part of the Kerguelen-Gaussberg Plateau 
(southern Indian Ocean). This archipelago is the subaerial continuation of Kerguelen hotspot 
volcanism for the last 30 Ma ( |Yang et al. ,"T998l |Weis" et al., 1998i [Nicolaysen et al., 2000D . 



The lava flows form the tabular reliefs (400 to 900 m high) observed today after glacial ero- 
sion, and represent more than 85% of the archipelago surface ( Giret, 1990D . The rest of the 



archipelago is composed of intrusions (gabbro, granite, and syenite) issued from Kerguelen 
plume melts QWeis &: Giret, 19*94 1 and quaternary glacial sediments. We studied paleomag- 



netic cores collected along three vertical sections of Kerguelen basalt at Mont des Ruches, 
Mont des Tempetes, and Mont Rabouillere sections (Fig.Q). For the reasons developed in Ple- 
nier et al. ()2()()2j] . the palaeomagnetic sections do not correspond exactly to the previously- 
dated sections ( Yang et al., 1998 Nicolaysen et al., 2000||Doucet et al., 2002 ), but they ap- 



pear to be correlated. Usually, we drilled seven samples in each successive lava flow using a 
gas-powered drill and oriented them with both solar sightings and magnetic compass with 
a clinometer. We took care to sample the bottom part of the least altered flows, and as far 
away as possible from intrusions. 
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3 ROCK MAGNETISM AND SAMPLE SELECTION 

For field intensities comparable to those of tlie Eartli's magnetic field (few tens of /iT), 
there is a proportionality between the Thermo-Remanent Magnetization (TRM) intensity 
measured at 20°C and the strength of the ancient magnetic field present during cooling 
through the blocking temperatures for almost all natural rocks. Thus, for some particular 
rocks cooling in the geomagnetic field during their formation, it is possible to estimate the 
palaeomagnetic field strength recorded by comparing their Natural Remanent Magnetization 
(NRM) with an artificial TRM acquired in the laboratory under a known ambient field. 
However, the coefficient of proportionality depends on grain size, shape distribution, and 
blocking temperatures as well as on the amount and type of ferromagnetic material the rock 
contains. In addition, this coefficient may have changed since the formation of the rock or 
during heating in the laboratory. For this last reason, a procedure using numerous successive 
heatings with increasing temperature steps has been developed ( |Thellier Thellier, 1959 1 in 



order to limit the field strength estimates to the temperature range preceding the irreversible 
magnetic and/or chemical changes in the ferromagnetic minerals. The strict conditions to 
be respected by the samples for correct palaeointensity determination are the following: 

(i) The Characteristic Remanent Magnetization (ChRM) recorded by the studied speci- 
men has to be a TRM, acquired at a known epoch in the geomagnetic field. 

(ii) The ChRM should not be disturbed by significant secondary magnetizations. 

(iii) The physical, chemical and crystallographical properties of the magnetic minerals 
must not have changed since the initial TRM was acquired nor changed during the successive 
heatings imposed by the experimental method. 

(iv) The independence and additivity laws of partial-TRM (pTRM) enunciated by Thel- 
lier ()1938|) have to be satisfied. That is, the total TRM must be equivalent to a sum of 
pTRMs, each associated with its own blocking temperature interval and not dependent on 
the remanence carried in every other intervals. This generally means that the magnetic 
carriers have to be single domain (SD) or in favorable cases pseudo-single domain (PSD) 
grains. 
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It is obvious that numerous samples can not fulfill these conditions, thus extensive prelimi- 
nary studies are necessary to avoid unnecessary work. 



3.1 Viscosity indices and demagnetizations 

We shown recently ( |Plenier et al., 2002D by means of a positive reversal test that the ChRM 
measured from Kerguelen lava are, generally speaking, not disturbed by unremoved sec- 
ondary components and then that these ChRMs are primary TRMs. In order to assess the 
importance of the secondary component carried by each sample, we analyzed the results from 
demagnetization experiments performed previously on sister specimens dPlenier et al., 2002 1 . 



First, we rejected samples for which the angle between the NRM and the ChRM was greater 
than 15° and those contaminated by a resistant secondary component (e.g., unremoved be- 
yond 20 mT alternating fields treatment or 300° C when samples are demagnetized by heating 
in zero field). Likewise, we kept only those flows for which the NRMs of all samples were well 
grouped. 60% of the flows fulfilled these conditions. For these flows, only the samples with a 
demagnetization curve as undisturbed as possible and presenting unblocking temperatures 
as high as possible, or median-demagnetization field of at least 20 mT, were retained. 
Second, we estimated the capacity of the specimens to gain a secondary Viscous Remanent 
Magnetization (VRM) by measuring them first after two weeks in the ambient magnetic field 
oriented along the z core axis and again, after two-week storage in a zero field. This enables 
determination of the viscosity index (Thellier fc Thellier, 1944) for each sample, which are 



reported in Table [TJ The viscosity index corresponds to about 25% of the VRM acquired in 
situ since the last reversal 780 ky ago ( |Prevot, 1981D . We choose the arbitrarily thresholds 
for the viscosity index of 10% as an upper limit for specimen from intermediate polarity flow 



and 5% for the others (Prevot et al., 1985). 



3.2 Susceptibility at room temperature and k-T curves 

To be sure of the thermal stability of the samples during successive heating in the labo- 
ratory, we used the low- field magnetic susceptibility (kg) measured at room temperature 
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after each thermal demagnetization step performed in air for the sister specimen previously 
studied in the palaeodirection determination dPlenier et al., 2002 1 . A favorable sample for 
palaeointensity determination should have a relatively constant kg value during most of 
the demagnetization procedure. However, this is not a sufficient criterion because thermally 
unstable samples may nonetheless display low variations in the susceptibility measured at 
room temperature. Thus to complete this approach, we measured continuously the low-field 
susceptibility of one sample from each fiow usually during two successive heating-cooling 
cycles under vacuum, the first up to 350°C and the second up to the Curie temperature 
(Tc). Fig. 121 presents two representative k-T curves (susceptibility as a function of tempera- 
ture) encountered during this study. The first case (Fig. |2t) illustrates the irreversible and 
complex thermomagnetic behaviour observed for almost 60% of the samples. The magnetic 
carriers are interpreted as original titanomagnetite associated with titanomaghemite, a prod- 
uct of their low temperature oxidation. We rejected the fiows yielding this behaviour because 
of their thermal instability. The second case (Fig. Eb) illustrates the reversible behaviour 
observed for the rest of the samples. The magnetic carriers are low-Ti titanomagnetites 
probably produced by high temperature oxyexsolution of the original titanomagnetites. We 
considered fiows presenting this second reversible behaviour as suitable for palaeointensity 
determination experiments. 

In order to complement this thermomagnetic investigation, we observed thin sections from 
each k-T curve type using an oil immersion objective. Fig. IHK and Eb show photomicro- 
graphs in natural light of sample 269b (fiow Tou2), which displayed irreversible behaviour. 
We saw an isotropic phase sometimes associated with a pleochroic ilmenite, as illustrated 
here. Because titanomagnetite and titanomaghemite are difficult to distinguish under the 
microscope, it is hard to conclude about the nature of the isotropic phase. These two miner- 
als are certainly both present, but the existence of cracks almost omnipresent in this phase 
suggest a larger amount of titanomaghemite. This interpretation agrees well with the ir- 
reversible k-T curve. Fig. Et and Eli show two minerals from a thermally stable sample 
(234c, fiow Teml6). They illustrate two different advanced stages of deuteric oxidation with 
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ilmenite or titanohematite lamellae exsolved from a residual titanomagnetite almost entirely 
altered. Again this observation confirms the interpretation of the k-T curves. 



3.3 Pilot analysis 

We kept only flows for which at least three samples from different cores fit all the previ- 
ously enumerated criteria for the palaeointensity experiments (32 out of the initial 57). We 
performed a pilot analysis with one or two samples from each of these flows in order to iden- 
tify the most favorable flows for reliable determinations and to define the more appropriate 
demagnetization steps for the next series. After applying a posteriori criteria which will be 
presented later, 5 flows from the Mont des Ruches section, 9 from the Mont des Tempetes 
section and 1 from the Mont de la Rabouillere section seemed suitable. Because a sample 
which did not pass all the a priori selection criteria may sometimes furnish a reliable deter- 



mination ( Coe, 1967b Perrin, 1998D , they should be regarded as serving only for selecting 



the most appropriate flows for palaeointensity determination. For this reason, the remaining 
samples from the suitable flows underscored by the pilot analysis have been incorporated in 
the following two determination series (68 specimens), including the samples which did not 
pass all the a priori criteria first. 

4 PALAEOINTENSITY EXPERIMENTS 
4.1 Experimental procedure 

We used the method of Thellier & Thellier ()1959|) in its classical form to estimate the 
palaeointensity of the geomagnetic field. We also followed a sliding pTRM checks pro- 
cedure dPrevot et al., 19*85 ) every two demagnetization steps in order to verify that the 



pTRM capacity remains unaltered when the heating temperature is progressively increased. 
Because of the nonlinearity of acquisition of TRM in multidomain grains ( |Levi, 19T7 ), ap- 



plying a laboratory field far from the recorded one may lead to underestimate the ancient 
field ( |Coe, 1967bt [lanaka fc Kono, 1984D . Thus, to improve the quality of data, we took 
care to apply, along the vertical axis of the samples, a constant field of 50 /xT (known with a 
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precision of 0.1 fiT), which corresponds approximately to the strength of the mean Ohgocene 
field for Kerguelen latitude. We also carried out the heating-cooling cycles under a vacuum 
better than 10~^ mbar to limit possible oxidation during experiments. We performed demag- 
netizations up to 580°C, with 13 steps ranging from 50 to 10°C, using a home made furnace 
(temperature reproduced within 2°C) in the palaeomagnetic laboratory of the University of 
Montpellier. Before the treatment, and after each heating-cooling cycle, we measured the 
remanence with a JR5-A spinner magnetometer. To ensure the reproducibility of the pro- 
cedure at the same demagnetization step (pair of heating-cooling cycles and pTRM check) , 
we always kept the specimens in the same place in the oven. 

4.2 Preliminary selection of palaeointensity data 

Because the interpretation of palaeointensity experiments is subjective, as many other data 
interpretations, a posteriori criteria help to ensure the technical quality and objectivity of 
the results and their comparison from study to study. It is noteworthy to point out that 
except the number of data used for each individual determination, which has been fixed 
to four consecutive points, the other a posteriori criteria were not used directly to reject a 
determination. When only one or two a posteriori criteria failed, we kept the determination 
choosing the nearest temperature interval for which the corresponding criteria stay as close 
as possible from the exclusion bounds fixed. The data with more than three independent 
criteria unfulfilled have been considered as unable to furnish a reliable palaeointensity esti- 
mate and thus have been excluded. 

(i) / criterion: A classical way of representing palaeointensity data is to use the Aral 
diagram dArai, 1963t [Nagata et al., 1963D in which the NRM^ remaining at each tempera- 
ture step T is plotted against the pTRM^ acquired in the laboratory field from T to room 
temperature. The slope of the least squares fit line computed from the linear part of the 
NRM-pTRM diagram gives an estimate of the palaeofield strength. Thus criteria have to 
be defined to constrain the determination of this best fit line and to quantify its technical 
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quality. 

We already fixed the minimum number of successive points used for the determinations to 
four. The other criteria used is the NRM fraction / given by the ratio of the NRM lost over 
the selected temperature interval to the total NRM ( |Coe et al., 1978 ). We consider that a 
value of 0.3 is the minimum NRM fraction for acceptable determination. Note that a quality 
factor q ( |Coe et al., 1978 Prevot et al., 1985) ) less than 5 can also indicate a determination 
of the palaeofield strength of poor quality. However, because q is not independent of /, this 
lead to reject the same samples in this study. Thus we did not use it as a rejection criterion. 



(ii) Control on vector endpoint diagrams (Zijderveld's plot): Because a straight line in 
the Aral plot can involve more than one component, we checked their existence on the 
Zijderveld's projection of the NRM demagnetization computed from the palaeointensity ex- 
periments. To complete this qualitative approach, we quantified the dispersion of the points 
regarding to the best fit line by the maximum angular deviation (MAD) ( |Kirschvink, 1980 1 
and chose a maximum value of 10° for this criterion. 

Another quality criterion is given by the angle a between the best fit line (anchored to the 
center of mass) and the vector average (anchored to the origin) of the selected data in the 
Zijderveld plot. If the data correspond to the primary remanence, a should be inferior to 
say 10°, in the opposite, they are certainly biased by a spurious unknown component. 

(iii) z ratio: The heating remanent magnetizations (HRM) ( |Calvo et al., 2002) ) acquired 
during the heating procedure under the application of a constant weak field in palaeointen- 
sity experiments, lead to erroneous data. Hopefully, the direction along which the laboratory 
field was applied corresponds to the axis of our cylindrical samples (Z axis). Therefore, in 
a Zijderveld plot in sample coordinates, the acquisition of HRM appears as a progressive 
deviation of the demagnetization curve in the vertical plane towards the vertical axis direc- 
tion. This is illustrated in Fig. |3 in which we compare the Zijderveld's plot inferred from 
the palaeointensity experiment for a rejected sample with the demagnetization curve of its 
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sister specimen obtained during the palaeodirection analysis. 

An "HRM check" is allowed using the ratio ( Goguitchaichvili et al., 1999a ): 

where HRMt and NRMt are the HRM created in the sample and the NRM left in the 
ChRM direction, at a given temperature T, respectively. Similar ratios exist (e.g. R or 
R' ( Coe et al., 1984| )) but the z ratio allows to monitor the evolution of alterations during 
treatment. Calculating this ratio supposes that one knows the ChRM direction, which re- 
quires demagnetization of a sister specimen before the palaeointensity experiment. 
To help the interpretation, we defined the upper limit of the accepted temperature interval 
as the one preceding a z >20%. However, we extended the interval over this limit when the 
ratio NRM left/TRM gained did not change in the new interval. In this case some magneto- 
chemical transformations are present but do not change the estimate. A limitation of the 
z ratio is that it depends on the angle between the ChRM and the applied field direction 
(direction of the HRM). The more this angle tends to 90°, the bigger is the deviation. Un- 
fortunately, the tray we used to place the specimens in the oven did not permit us to orient 
their ChRM at 90° to the applied field. Hence, the z we calculated is a minimum value, 
thus these steps had to be rejected. Another limitation of this criterion is that some unex- 
pected fluctuations may appear at high temperature when the NRM left is too small. For 
this reason, we completed this approach by observing the evolution of the demagnetization 
on an equal area projection. In Fig. IH the z-T curve of a rejected sample 163E (Tem06) is 
compared to the one obtained from an accepted sample (213E, flow Teml3). 

(iv) Difference ratio (DRAT): We performed sliding pTRM checks in order to estimate 
the temperature at which alteration of the magnetic minerals begins. Commonly, many 
authors consider that repeated pTRM acquisition at the same temperature steps should 
agree to within 15% dGoguitchaichvili et al., 1999bD . However, pTRMs acquired within low- 
temperature interval, are rather small and thus a reproducibility within a given percentage is 
difficult to obtain. Therefore, we used the difference ratio (DRAT) dSelkin fc Tauxe, 2000 ), 
which normalizes the maximum difference between repeated pTRM steps performed on the 
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temperature interval selected for the determination by the length of the corresponding NRM- 
pTRM segment. We fixed a maximum acceptable value of 10% for this criterion. However, 
this ratio reveals possible alterations on the estimate interval only, whereas physical and 
chemical alterations of the heated samples can appear as soon as the first demagnetiza- 
tion steps ( Kosterov fc Prevot, 1998) ). Therefore, we also monitored the lower temperature 



pTRM checks by determining the DRAT from ambient to maximum temperature of the 
interval used to estimate the palaeointensity. It is noteworthy that the sample 163E (Fig.^, 
rejected because of its z ratio, has a DRAT of only 3%. Also five samples which give an esti- 
mate in Table m possess a DRAT >10%, thus these two alteration criteria are complementary. 

(v) High temperature pTRM-tail test: Unfortunately, 56% of the samples present a more 
or less pronounced curvature in their NRM-TRM diagrams and some specimens provide two 
acceptable estimates regarding the selection criteria. Because independence and additivity 



laws of pTRM ( Thellier, 1938 1 are violated for MD but maybe also for PSD grains, the 
determination of the domain structure could be a relevant tool to choose between the two 
acceptable palaeointensity estimates. 

In case of thermally stable samples, we used the thermomagnetic criterion that was first 
introduced by Bol'shakov & Shcherbakova ()1979p to get insight on the domain structure of 
ferrimagnetics. We calculated the parameter Aht which is illustrated in Fig. El and defined 
by: 

= ta.lHT|pTRM(T. T.)| 

^' ^' pTRM(Ti,T2) ^ ' 

where pTRM(Ti,T2) is the pTRM measured at room temperature gained between Ti and 
T2 (Ti < T2) during cooling from the Curie temperature (Tc) in a lOO^uT field, and the 
tailHT[pTRM(Ti, T2)] is the part of this pTRM not demagnetized when the sample is sub- 
sequently heated again to Ti and cooled down is zero field. According to the criteria defined 
by Shcherbakova et al. ()2()()()|1 for Aht(Ti, T2) <4%, the remanent carrier are predominantly 
SD grains, for 4% <Aht(Ti,T2) <(15-20)%, they behave as pseudo-single domain (PSD) 
grains, and for Aht(Ti,T2) >20% they are predominantly MD grains. 
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We measured the coefficients Aht(Ti,T2) at increasing temperature intervals for one ther- 
mally stable sample per flow in order to determine the best temperature interval to estimate 
the paleointensity. Because the heatings were performed in air, we repeated the determina- 
tion of the coefficient Aht(300, Troom) two times, at the beginning and at the end of the 
treatment, in order to monitor alterations appearing during the successive heatings. Unfor- 
tunately, the vibrating thermal magnetometer (VTM) we used broke before we could mea- 
sure all the representative samples of each selected flow. The results available are reported 
in Table |21 Even though all the flows have not been studied, the global trend observed 
confirms that with increasing temperature, the magnetic carrier behaviour is more SD- 
like ( Carvallo et al., 2003[ [Shcherbakova et al., 2000 ). Thus, when two coexisting palaeoin- 



tensity estimates remain possible (5 samples), the high temperature pTRM-tail test favors 
acceptance of the higher temperature interval determination, as for sample 187E shown in 
Fig. El 

(vi) Low temperature pTRM-tail test: According to Fabian ()200ip . the pTRM tail is not a 
direct measure for a sample's tendency to yield a curved Aral plot. The common concave-up 
shape of MD Aral diagrams is rather well accounted for by the low temperature tail (Low-T 
tail) of each pTRM segment acquired during cooling from Tc (pTRMlt(Ti, T2)). This Low- 
T tail is due to the grains with unblocking temperature less than blocking temperature. In 
order to compare the results obtained from the two tail tests and validate the choice of the 
temperature intervals for the interpretation, we decided to evaluate the Low-T tail of some 
characteristic samples used during the palaeo intensity experiment. We define the coefficient: 

tailLT[pTRM(Ti,T2)] 
- pTRM(Ti,T2) ^ 

where pTRM(Ti, T2) is the pTRM measured at room temperature gained between Ti and T2 
(Ti < T2) in a 50/iT field during cooling from the Curie temperature and tailLT[pTRM(Ti, T2)] 
is the part of this pTRM removed when the sample is subsequently heated to T2 and cooled 
down in zero field (Fig. EJ. This coefficient is analogous to the Aht coefficient defined by 
Shcherbakov et al. ()200ip but the relation between the magnitude of the two tails is not clear. 
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We can expect from Dunlop & Ozdemir (|2U(Jip that Alt is equal to Aht because the un- 
blocking temperature distributions are almost symmetric for a given blocking temperature, 
but this may be quite different for hundred- degree temperature intervals. In the absence of 
a quantitative study, only the evolution of the Alt ratio for increasing temperature intervals 
will be considered. 

The Alt coefficients for the (450,350) and (550,450) temperature intervals are reported in 
Table El The main result is that Alt(550,450) is always smaller than Alt(450, 350); thus 
the Low-T tail test again favors the higher temperature interval in case of two acceptable 
determinations. In order to estimate the alteration which may occur during the required pre- 
liminary heating up to T^, the untreated samples 187D and 201C have been incorporated in 
the experiment. Comparing their results with the data from the corresponding sister samples 
187E and 201E, we can observe that, even though these samples appear thermally stable 
(reversible k-T curves), Alt(450,350) is greater for the sample used in the palaeointensity 
estimate and on the contrary Alt(550,450) is greater for the "virgin" sample. Nonetheless, 
for the five samples that allowed two acceptable interpretations, we choose the higher tem- 
perature interval, which yields determination more coherent with the other samples from 
the same flow. 

4.3 Palaeointensity results 

Out of the 402 samples from the Mont des Ruches, Mont des Tempetes, and Mont Rabouillere 
sections, 102 have been chosen using a priori criteria. After selection following a posteriori 
criteria, 49 samples furnish a technically acceptable determination (Tabled}. 
The quality factor q varies from 1.6 to 54.3 with a mean value of 12.1 for estimates made 
on a mean of 7 points. However, some a angles are relatively large, which could be due to 
an overprint of viscous origin, even though the samples are not strongly viscous, or to a 
secondary component acquired during the treatment, as shown by the angle (3 between the 
ChRMs of the sample and its sister specimen. Although secondary HRM exists, we consider 
its effect negligible on the determination, considering the z criterion, the linearity of the Aral 
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diagram and the evolution of the demagnetization in an equal area projection. Therefore, 
as evidenced by the MAD, which is the most widely satisfied a posteriori criterion, we are 
confident about these interpretations. Moreover, as is commonly done in palaeointensity 
studies, we distinguished two types of determination: class A for the samples which passed 
successfully all the criteria and class B for those in which one to three criteria were unfulfilled. 
However, as observed with flows Rabl2 or Tem6, there is no single relation between the 
field estimate and its classification. Thus, the interpretations made with a few unsatisfied 
a posteriori criteria can be also considered as reliable. Fig. [7| shows some characteristic 
examples of palaeointensity determinations. 

We present, in Table EJ the mean palaeofield strength with its associated standard deviation 
for each individual flow. Three consecutive Mont des Ruches flows (Rucl4, 15, and 16) yield 
a well defined palaeointensity in the range 35-40 fiT. This value is comparable with the 
only reliable flow from the Mont Rabouillere section (Rabl2) and with the flows Teml3 and 
Teml5 from the Mont des Tempetes. For this last section we also observe three significantly 
higher values around 50-60 /iT (Teml,17 and 18) and two well defined lower field strengths 
of only 20.1 and 25. Oil. 7 /iT (Tem6 and 10). Because, the within flow standard deviation 
represent at the maximum 15% of the palaeointensity value (flow Teml3), we are quite 
confident about the reliability of these estimations. 

We performed in a previous study (|Pl enier et al., 2002] ), a quantitative bootstrap test for 
a common mean on successive magnetic directions in order to identify non-independent 
records. We concluded in particular that the Rucl4 and Rucl5 flows may represent two 
contemporaneous records of the palaeofield (TableU}. The palaeointensity estimates obtained 
for these two successive flows overlap within their uncertainties and thus strengthen our 
former interpretation. However, because the palaeomagnetic field may remain constant for 
a relatively long time ( |Love, 2000| ) we will nevertheless consider in the present study each 
of the 12 reliable palaeointensity estimates as distinct records. 

Concerning the samples lOlF, lllF, 201E, 160E, 161F and 229C, for which pTRM-tail tests 
were performed, we decided to keep the lower rather than the higher temperature interval 
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estimates because the latter was sometimes impossible to define owing to alteration of the 
pTRM blocking temperature spectra (lOlF, 160E, 201E and 229C) or meaningless compared 
to the other determinations from the same flow (lllF, 161F). 

5 DISCUSSION 

5.1 Toward an improvement of palaeointensity determination 

Because of the strict conditions imposed by the commonly used palaeointensity determina- 
tion method of Thellier & Thellier ()1959p . numerous checks are needed to ensure reliable 
estimates. For this reason, palaeointensity experiments are time consuming and at the end 
many samples do not furnish reliable determinations. Consequently, the number of reliable 
estimates is very low and studies dealing with palaeointensity determinations are relatively 
rare. Then, effective selection followed by careful palaeointensity determination method are 
desirable to increase the productivity of the experiments and largely implement the number 
of palaeointensity estimates in the databases. 

However, no ideal selection criteria exist to select suitable specimens and flows for palaeoin- 
tensity experiments. For example, Sherbakov et al. ()2()()1|1 proposed to use a pTRM-tail 
check as a systematic procedure before the experiments to discriminate the samples with 
MD grains. The few determinations of the low and high temperature pTRM-tail performed 
in the present study are not convincing since many samples yielding a paleointensity estimate 
with good technical quality (Table |^ will not have been kept if we followed Shcherbakov et 
al.'s (l2UUip recommendation. Thus, we do not recommend to use the thermomagnetic crite- 
rion for selection but rather in a case to case basis to define the best temperature interval 
for thermally stable specimens exhibiting two slopes in the NRM/TRM diagram. For the 
selection based on k-T curve shapes, the temperature intervals used for the estimates (c.f. 
Table 13)) can change at the flow scale (flow Teml3 or Teml5) without important deviation 
in the palaeofield strength. Therefore, the thermal behaviour of the samples varies within a 
flow and this criterion applied to a supposedly representative sample is not totally effective. 
A study of each successful sample, at the end of the palaeofield strength determination pro- 
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cedure, up to the maximum temperature reached for the palaeointensity estimate would be 
more interesting. Whatever the numerous a priori precautions we took in this study, only 
48% of the selected samples gave a reliable estimate. Thus, in the absence of relevant a priori 
criteria, it is preferable to reduce the preliminary experiments to the determination of the 
viscosity index, which is an already rapid and reliable way to discriminate the specimens 
disturbed by a viscous secondary component, followed by a directional analysis, needed for 
careful palaeointensity estimates and giving enough information to discard possibly prob- 
lematic flows because of large NRM scatter, poor technical quality of demagnetization curves 
and so on. 

As regard to the determination procedure itself, some modifications of the Thellier's method 
have been proposed recently ( |Calvo et al., 2002 Riisager fc Riisager, 2001D . They include 
an additional heating at each demagnetization step in order to control the HRM creation 
and/or a "pTRM tail check". However, the lack of suitable selection criteria lead to more 
failure of the palaeointensity experiment and these more time consuming procedures are 
consequently not helpful for systematic studies. Orienting the samples so that their ChRM 
is perpendicular to the direction of the applied field is less time consuming and sufficient 
to check for possible HRM creation. Likewise, at the end of the treatment, complementary 
experiments like k-T curve shape determination, extended Thellier's method ( |Fabian, 2001| ) 
for thermally stable samples or "pTRM tail checks" can be performed to ensure the qual- 
ity and sharpen the reliable estimates, but such supplementary processes must concern the 
acceptable samples only. Thus, good determinations do not necessarily need supplementary 
heatings nor burdensome procedures to be performed; this should be limited only to the 
problematic samples. Reducing the selection and placing the control procedures at the end 
of the treatment is thus a way to perform reliable palaeointensity estimates more quickly, 
without lowering the quality of the determinations. 
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5.2 Comparison with previous palaeointensity results 

Five Oligocene flows from tlie He Haute section (49.4°S, 69.9°E) liave been processed for 
palaeointensity study ( |Derder et al., 1990 ). Of the 24 samples analyzed, 11 yielded a deter- 



mination, but only one flow provided more than three estimates of the palaeofield strength. 
Moreover, the uncertainty of the mean palaeointensity is equal to 36% and consequently, 
these data cannot be considered as reliable. The 12 distinct mean flow estimates presented in 
this paper are thus the only reliable determinations available for the Kerguelen Archipelago. 
However, data well distributed on the Earth's surface are needed to provide a correct idea 
of the palaeomagnetic field behaviour. It is then important to pool the Oligocene Kerguelen 
results with other reliable determinations already achieved. For this aim, we used the up- 
dated lAGA 2002 palaeointensity database available at the following address: 

ftp : //ftp . dst u. univ-montp2 . fr / pub / paleoint db / 
This database presents estimates obtained with different quality determinations. However, 
inclusion of lower quality data leads to higher average values of the geomagnetic field 
( Juarez fc Tauxe, 2000| [Goguitchaichvili et al., 1999b ), thus a selection using identical cri- 



teria to those used in this study is needed before we can combine them meaningfully with 
the Kerguelen results. For this reason we discarded estimates based on intermediate polar- 
ity flows, obtained with other methods than the Thellier & Thellier's f)1959|) associated 
with pTRM checks, presenting less than 3 individual samples from the same unit (except 
for basaltic submarine glasses) and a standard deviation of the mean >20%. Application of 
these drastic selection criteria led to the rejection of 660 determinations out of the 910 esti- 
mates contained in the database between 0.3 and 50 Ma (Fig. |H^). Then we sharpened the 
selection to a time window between 20 and 40 Ma, which resulted in 11 suitable estimates 
among the 74 initial extended Oligocene records of the data base (Fig. |HK)- 
In order to compare the Kerguelen determinations with the selected estimates issued from 
various locations and recording different time, we calculated the virtual dipole moment 
(VDM) given by: 



18 G. Plenier, P. Camps, R.S. Coe and M. Perrin 

VDM = ^Fancient(l + 3 COS^ 9)-'/' (4) 

which corresponds to the moment of a dipole field producing the estimated palaeointensity 
Fancient at the magnetic colatitude 6. R is the Earth radius. Results of the calculations are 
reported in Table 0] and the 12 suitable determinations are plotted in Fig. |H| 
Even though we performed careful estimates in order to avoid any possible MD effects in 
the interpretations, which could lead to systematic higher determinations, the arithmetic 
mean obtained with the Kerguelen data, 6.15±2.1 10"^"^ Am^, is definitely higher than the 
other Oligocene estimates already achieved ( Goguitchaichvili et al., 200H [Riisager, 1999 



Juarez et al., 1998 ). Including Kerguelen determinations leads to an increase of the Oligocene 
mean VDM from 4.1±0.5 to 5.4±2.3 10^^ Am^. The new Oligocene mean VDM is closer to 
the mean VDM estimated for the 0.3 and 5 Ma interval Puarez fc Tauxe, 2000 1 and higher 
than the mean VDM defined between 0.3 and 300 Ma (Fig.|Hl). Thus these selected Oligocene 
estimates favor a relatively stable field between 0.3 and 40 Ma and reinforce the idea of an 
exceptionally high recent geomagnetic field strength. However, the lack of good palaeoin- 
tensity data didn't allow us to go further in our interpretation. As illustrated in Fig. |Ht, 
the Oligocene data do not uniformly cover the Earth's surface even though the Kerguelen 
results equalize the number of estimates available for each hemisphere. More reliable deter- 
minations are needed to better constrain the evolution of the palaeofield with geological time. 
It is important to note that the Kerguelen results represent half of the reliable palaeointen- 
sity estimates available for the Oligocene, which could lead to a bias towards higher values. 
Moreover, only 3 extended Oligocene determinations are from normal polarity flows. Thus 
even if the 12 new estimates are supposed to be independent, the Ohgocene palaeomagnetic 
field may not be sufficiently time-averaged. 



6 CONCLUSION 

We carried out a palaeointensity study on three Oligocene (28-30 Ma) volcanic sections 
from the Kerguelen Archipelago (southern Indian Ocean). It complements the palaeosecular 



Absolute palaeointensity from the Kerguelen Archipelago 19 
variation directional study already achieved on the same sections: Mont des Ruches, Mont 
des Tempetes, and Mont Rabouillere dPlenier et al., 2002 ). Out of the 57 studied units, 
we considered 32 suitable for palaeofield strength determination experiments when at least 
three of their samples pass the following selection criteria: angle between NRM and ChRM 

< 15° (or secondary component quickly demagnetized), NRM dispersion not too large, vis- 
cosity index <5 %, stable susceptibility at room temperature after each demagnetization 
step of a sister sample and reversible susceptibility in temperature curves. After preliminary 
experiments, we studied all samples from 12 favorable flows. We considered only the deter- 
minations made with at least 4 successive steps, a z ratio ( Goguitchaichvili et al., 1999a I 

< 20%, a DRAT dSelkin fc Tauxe, 2000D < 10%, f factors (|Coe et al., 1978|) > 0.3, a MAD 



< 10° ( Kirschvink, 1980 ) and an angle a between the best fit line and the vector average 

< 10° to be of good technical quality. In cases of two acceptable estimates, pTRM tail 
tests were used to choose the temperature interval corresponding to the more SD like be- 
haviour. This careful interpretation of the data leads to 49 reliable estimates. The VDMs 
calculated for the 12 flows vary from 2.78 to 9.47 with an arithmetic mean value of 6.15±2.1 
10^2 Am^. This study gives the first reliable palaeointensity estimates from the Kerguelen 
Archipelago and significantly increases the number of Oligocene data of comparable qual- 
ity. The new Oligocene mean VDM calculated, 5.4±2.3 lO^^Am^, is very close to the 0.3-5 
Ma value of Juarez & Tauxe ()2000p (5.5±2.4 lO^^Am^), and suggests little evolution of the 
geomagnetic field between 0.3 and at least 40 Ma. Thus the present moment of the field (8 
10^2 Am^) is confirmed to be exceptionally high. However, the lack of reliable data limits the 
interpretation and we discuss practical solutions to speed the selection of suitable samples 
for careful determination procedures in order to facihtate systematic studies and increase 
existing palaeointensity databases with high-quality estimates. 
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Figure 1. Location of the studied sections: a) Mont des Ruches : 18 flows (48°52'18"S, 68°54'48"E), 
b) Mont des Tempetcs : 20 flows (48°52'50"S, 69°06'37"E), c) Mont Rabouillere : 19 flows 
(49°05'25"S, 69°26'25"E). 
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Figure 2. Low-field susceptibility versus temperature curves: continuous measurements at temper- 
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Table 1. Cleaned average directions of magnetization of lava flows retained for paleointensity exper- 
iments from Plenier et al. (2002). 



Flow Chron n/N Inc Dec 095 k Plat Plong i/% Alt 



Mont des Ruches section (48.87°S, 68.91°E) 



Rucl6 


C9r 


4/4 


79.9 


183.3 


6.5 


201.6 


-68.4 


65. 


,9 


2, 


.48 


180 


Rucl5(i) 


C9r 


7/7 


74.5 


159.2 


3.9 


244.8 


-73.1 


105. 


,3 


4, 


.60 


160 


Rucl4(i) 


C9r 


7/7 


76.4 


149.5 


4.5 


183.8 


-67.7 


104. 


,5 


5, 


.80 


150 






Mont 


des Tempetes section (48. 


88°S, 69 


.11°E) 












Teml8 


C9r 


8/8 


60.7 


177.0 


3.5 


248.7 


-82.5 


231. 


,6 


2, 


.73 


172 


Teml7 


C9r 


7/7 


62.0 


186.0 


4.0 


224.8 


-83.0 


287. 


,7 


3, 


.36 


160 


Teml6 


C9r 


5/7 


64.2 


178.1 


4.2 


328.6 


-86.8 


224. 


,6 


5, 


.57 


152 


Teml5 


C9r 


4/4 


64.5 


188.1 


2.8 1109.8 


-84.0 


317. 


,3 


4, 


.04 


145 


Teml3 


C9r 


7/7 


55.5 


180.5 


5.0 


146.0 


-77.2 


250. 


,9 


2, 


.99 


130 


TemlO 


C9r 


7/7 


60.6 


188.0 


7.0 


75.6 


-80.8 


289. 


,7 


2, 


.56 


105 


Tem6 


C9r 


5/5 


77.1 


221.2 


2.4 


979.6 


-63.0 


31. 


,9 


4, 


.94 


75 


Teml 


ClOn.ln 


7/7 


-67.3 


356.6 


2.8 


464.7 


87.5 


309, 


,2 


5, 


.87 


5 






Mont Rabouillere section (49.09°S, 69.44°E) 












Rabl2 


ClOr 


7/7 


68.2 


140.2 


4.2 


207.0 


-64.8 


138. 


,7 


1, 


.80 


150 



indicates flows which have been grouped together in Plenier et aL's (2002) directional analysis. Flows are listed in stratigraphic 
order with the youngest on top, oldest on the bottom. Chron correspond to the polarity chrons inferred from Plenier et al.'s (2002) 
analysis. n/N is the number of samples analyzed/total number of samples collected. Inc and Dec are the mean inclination, positive 
downward, and the declination east of north, respectively. 095 is the 95% confidence envelope for the average direction, k is the 
precision parameter of Fisher distribution. Plat/Plong is the latitude/longitude of VGP position, respectively. u% is the geometric 
mean viscosity index j'l'hellier &: 'I'hellier, 1944} . Alt is the altitude of the flow in meter. 



Table 2. High-temperature pTRM-tail test 



Flow 


Sample 


SOO-Tt-ooiti 


400-300° C 


500-400° C 


550-500° C 


30n-T* 






Aht (B) 


Aht (B) 


Aht (B) 


Aht (B) 


Aht 


Rucl6 


117B 


17.0 


34.6 


n.d. 


n.d. 


n.d. 


Rucl5 


lllC 


48.8 (9.7) 


48.1 (11.5) 


30.6 (23.3) 


13.9 (55.5) 


50.7 


Rucl4 


lOlC 


20.8 (11.7) 


29.8 (7.5) 


13.3 (27.9) 


4.0 (52.9) 


17.3 


RuclO 


076B 


25.2 (40.6) 


33.1 (18.1) 


32.7 (12.9) 


6.8 (28.4) 


20.3 


Teml3 


216C 


14.4 (11.4) 


30.1 (9.9) 


9.5 (32.4) 


10.3 (46.3) 


21.1 


TemlO 


187B 


23.9 (10.4) 


25.5 (10.7) 


11.6 (31.8) 


5.0 (47.1) 


21.3 




189B 


29.8 (10.1) 


34.1 (7.3) 


21.7 (19.1) 


5.3 (63.5) 


26.4 


Tem6 


160B 


32.8 (14.9) 


17.4 (16.6) 


28.6 (20.1) 


11.5 (48.4) 


30.7 


Teml 


201B 


12.4 (17.5) 


15.5 (3.0) 


16.6 (21.3) 


3.7 (58.2) 


8.8 


RablO 


620D 


15.2 (17.1) 


35.0 (6.7) 


11.4 (34.8) 


7.6 (41.4) 


16.8 



Kht values are the relative intensities measured at room temperature of the high temperature 



pTRM-tafl expressed in percent A(Ti,T2) = tan[pTRM(Ti, T2)]/pTR,M(Ti, T2). B values shown 
in parentheses correspond to the percent of the total pTRM (e.g., pTRM;) each pTRM(Ti,T2) 
represents; nd means not determined. * is a pTRM-tail test (300-T,-oom) repeated at the end of the 
treatment to control the thermal stability of the sample. For Aht(Ti , T2) <4%, the remanent carriers 
are predominantly SD grains, for 4% <Aht(Ti,T2) <(15-20)%, they present a PSD behaviour, and 
for Aht(Ti,T2) >20% they are predominantly MD grains. 
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Figure 3. Reflected light microphotographs using an oil immersion objective. a),b): sample 269b 
(Tou2) in natural light shown in different orientations in order to illustrated the pleochroic ilmenite 
associated with a non pleochroic phase which could correspond to a titanomaghemite; c),d): two 
characteristic minerals from sample 234c (Teml6) in polarized light. The crossed nicols are at 90° 
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Figure 5. Continuous thermal demagnetization of pTRM(Ti,T2). The low-temperature pTRM- 
tail for PSD (A) or MD (B) grains corresponds to the part of pTRM(Ti,T2) removed at T2, while 
the high-temperature pTRM-tail corresponds to the part of this pTRM unremoved at Ti, (C) for 
MD and (D) for PSD grains. The low and high-temperature pTRM-tails are measured at room 
temperature. 
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Figure 6. Arai diagram for sample 187E (flow TemlO) illustrating a case for which two acceptable 
interpretations are possible. A) The low-temperature interval is rejected regarding the pTRM- 
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Figure 7. Characteristic examples of palaeointensity estimates from this study, a) and b) corre- 
spond to grade A estimates, c), d), e) and f) to grade B estimates. Black (White) circles are the 
step used (not used) for the palaeointensity estimate, the triangles represent the pTRM checks. 
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Table 3. Low-temperature pTRM-tail test 



Flow 




450-350°C 


550-450° C 






A.TT fB) 




Teml8 


251C 


15.7 (9.0) 


2.9 (64.6) 


Teml6 


237E 


8.2 (9.3) 


2.0 (74.1) 


Teml5 


229C 


12.3 (9.4) 


5.1 (59.8) 






t .o yo.^t J 




Teml3 


213E 


6.9 (7.6) 


2.3 (31.8) 




214D 


23.5 (7.7) 


3.9 (18.8) 




215F 


3.6 (8.1) 


3.4 (22.6) 




217D 


5.7 (6.5) 


2.9 (43.4) 


TemlO 


187E 


13.0 (7.0) 


2.3 (65.1) 




187D 


8.8 (7.2) 


2.4 (54.8) 




189E 


12.9 (5.1) 


2.3 (57.9) 


Tem6 


161F 


3.8 (11.7) 


1.5 (42.9) 


Teml 


201E 


18.3 (3.3) 


2.2 (72.0) 




201C 


10.1 (6.3) 


2.4 (73.9) 



Alt(Ti,T2) are the relative intensities measured at room 



temperature of the part of the pTRM(Ti,T2) removed after 
heating to T2 in zero field. B values shown in parentheses cor- 
respond to the percent of the total TRM (pTRM(580,r^t,om)) 
each pTRM(Ti,T2) represents. 
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Table 4. Palaeointensity determinations 



Teml8 



Teml7 



Tern 16 

Teml5 
Teml3 



TemlO 

Tem6 

Teml 

Rabl2 



Spl 


Grade 


FeibcrFe 


AT 


n 




f 




g 


q 


MAD 


a 


/? 


DRAT 


115C 


B 


34.0±1, 


.4 


20-400 


5 


0, 


,32 


0, 


,72 


5.5 


6.2 


(11.7) 


7.9 


0.7 


116C 


B 


34.2±2, 


.6 


150-400 


4 


(0, 


,22) 


0. 


,58 


1.6 


3.5 


7.9 


5.7 


0.4 


117F 


B 


QA 8^9 


9 


150-400 


5 


(0, 


,23) 


0, 


,62 


2.2 


2.2 


7.5 


3.0 


4.1 


118F 


B 


OO.OILl, 


Q 
.O 


20-400 


5 


0, 


,32 


0, 


,68 


6.2 


1.7 


(16.6) 
8.8 


14.8 


3.5 


108C 


B 




Q 


20-400 


5 


(0, 


,25) 


0, 


,73 


2.7 


4.2 


7.8 


0.4 


llOE 


B 


Q7 0_|_1 


1 

. 1 


350-540 


8 


0, 


,37 


0, 


,81 


10.2 


4.4 


(19.8) 


16.7 


4.1 


lllF 


B 


4^ n+i 


T 

. X 


20-400 


6 


(0, 


,23) 


0, 


,79 


4.9 


4.6 


11.2 


10.8 


1.7 


114D 


B 


QA 9+n 


Q 


150-480 


6 


0, 


,53 


0, 


,79 


15.9 


1.8 


(16.4) 


11.0 


5.2 


lOlF 


B 


QQ p,-L-\ 


1 

. 1 


20-480 


8 


0, 


,37 


0, 


,81 


10.8 


4.1 


(15.7) 


12.9 


2.7 


102D 


B 




n 


150-495 


7 


0, 


,39 


0, 


,80 


12.6 


2.5 


(12.2) 


10.5 


2.4 


103D 


B 


oo. Ozti, 


Q 
.O 


300-495 


6 


0, 


,33 


0, 


,75 


4.5 


2.2 


(12.1) 


5.9 


8.1 


245D 


A 


65.4±2, 


.1 


20-510 


9 


0, 


,44 


0. 


,86 


11.8 


4.2 


8.2 


2.6 


5.5 


246C 


B 


43.9±0, 


.8 


300-495 


6 


(0, 


,27) 


0. 


,78 


11.5 


2.7 


(18.2) 


10.1 


7.6 


248D 


B 


•^1 fi+1 


8 
.O 


150-500 


8 


(0, 


,29) 


0. 


,82 


7.2 


5.7 


(13.0) 


10.0 


4.2 


249D 


A 


r.o n_(_4 


.2 


20-400 


5 


0, 


,36 


0. 


,71 


3.2 


2.6 


7.3 


1.9 


0.1 


251C 


A 


54.8±1, 


.5 


150-525 


9 


0, 


,51 


0. 


,85 


15.8 


3.2 


4.9 


3.9 


4.1 


252D 


B 


fil 9+1 


fi 
.u 


20-495 


8 


0, 


,52 


0. 


,81 


16.0 


3.0 


(10.1) 


6.4 


(13.8) 


238D 


B 


/la /I4-9 

4o.4ztZ, 


.0 


150-495 


7 


(0, 


,28) 


0, 


,82 


4.5 


3.4 


(14.7) 


11.5 


5.7 


239D 


B 


58.8±z, 


-1 

.i 


150-480 


6 


(0, 


,25) 


0, 


,78 


5.4 


5.1 


(19.4) 


14.0 


(10.1) 


240E 


B 


3y.z±o, 


.9 


150-500 


8 


0, 


,40 


0. 


,79 


13.1 


5.1 


(11.1) 


5.7 


5.3 


241D 


B 


56.6±1, 


.1 


150-495 


7 


0, 


,36 


0. 


,82 


15.1 


3.3 


(12.5) 


7.2 


6.3 


243E 


B 


57.3±3, 


.5 


150-450 


5 


fo 


,23) 


0, 


,73 


2.7 


4.9 


(26.3) 


27.5 


6.6 


244E 


B 


50.3±0, 


.9 


1 50-51 


« 

O 


n 


35 


n 


84 


1 6 


5 


n 1 


5 9 


2.2 


232D 


A 


35.2±0, 


.8 


20-540 


11 


n 


54 





87 


90 8 


9 8 


8 


1 


4.2 




T3 
JD 


30.4±0, 


.4 




c 
o 


n 


,uy 


n 


84 

, 04 




9 


9 


1 Q 


(^ n 7'\ 




T3 
D 


24.6±0, 


.4 


OUU-D4U 





n 


. 1 6 


u. 


8/1 
,o4 


oO.o 


0.4 


i.o 


Q n 
o.U 


(iti.lj 


Zo / th 


A 

/\ 


31.6±0, 


.6 


O0U-04U 





n 


,04 


n 
u. 


,00 


18 1 
io. 1 


9 <^ 
Z.O 


1 

i.o 


A 

4.0 


9 

Z.O 




A 

JrV 


34.3±0, 


.8 


4sn ^sn 


e 
o 


n 


7fi 


n 


, t u 


9^ 


1 


1 S 


9 fi 


7 7 


99QP 




38.2±1, 


.7 


1 i^n zLQ^ 
iou-4yo 


7 


n 


'iC\ 
,ou 


n 


89 
,OZ 


^ A 

0.4 


7 9 
1 .Z 




^'^ A 

iO.4 


9 1 
Z. 1 




A 


37.1±0, 


.4 


480-580 




n 


78 


n 


71 


54.3 


1.6 


0.3 


5 1 


5 8 


213E 


A 


32.8±1, 


.0 


525-580 


5 


0, 


,41 


0. 


,72 


9.8 


2.2 


0.9 


2.8 


2.0 


214D 


A 


34.7±1, 


.7 


20-555 


12 


0, 


,48 


0. 


,85 


8.1 


3.5 


5.8 


6.5 


2.5 


215F 


A 


36.6±1, 


.5 


20-555 


12 


0, 


,58 


0. 


,85 


12.4 


5.1 


2.2 


4.6 


2.4 


216F 


A 


40.0±1, 


.9 


20-500 


9 


0, 


,52 


0. 


,85 


9.1 


7.8 


8.5 


4.0 


5.4 


217D 


A 


45.9±3, 


.3 


20-555 


12 


0, 


,74 


0. 


,85 


8.7 


1.9 


3.3 


2.7 


3.2 


218C 


A 


48.3±3, 


.1 


20-555 


12 


0, 


,50 


0. 


,82 


6.3 


3.5 


6.4 


7.3 


5.2 


187E 


A 


9zL 8-1-1 

Z^.OZIZ 1 , 


n 


480-540 


5 


0, 


,35 


0. 


,71 


6.3 


2.2 


5.4 


2.5 


5.5 


188E 


A 


27.2±1, 


.4 


510-555 


4 


0, 


,47 


0. 


,53 


4.8 


4.1 


4.7 


3.2 


5.1 


189E 


A 


23.0±0, 


.5 


525-570 


4 


0, 


,66 


0. 


,62 


17.5 


3.4 


1.1 


6.3 


6.1 


160E 


B 


21.4±0, 


.5 


20-500 


9 


0, 


,57 


0. 


,83 


21.4 


9.4 


(14.9) 


6.7 


6.0 


161F 


A 


21.2±0, 


.9 


150-525 


9 


0, 


,53 


0. 


,85 


10.6 


3.6 


6.0 


2.8 


4.0 


162E 


B 


17.8±0, 


.5 


20-525 


10 


0, 


,52 


0. 


,80 


16.7 (10.6) (14.0) 


6.1 


5.5 


201E 


B 


55.9±4, 


.9 


150-480 


7 


0, 


,31 


0. 


,82 


2.9 


6.1 


(16.9) 


13.2 


1.8 


203E 


B 


62.9±2, 


.6 


20-495 


8 


0, 


,42 


0. 


,84 


8.6 


3.8 


(11.0) 


9.8 


7.3 


205E 


B 


64.1±3, 


.2 


150-510 


8 


0, 


,35 


0. 


,82 


5.7 


4.7 


9.9 


6.6 


(13.1) 


628B 


A 


35.5±0, 


.6 


225-500 


7 


0, 


,53 


0. 


,77 


24.4 


3.3 


5.7 


4.3 


5.2 


629B 


A 


38.0±2, 


.0 


150-480 


7 


0, 


,54 


0. 


,75 


7.7 


2.1 


8.1 


5.2 


7.1 


630C 


A 


41.4±1, 


.1 


150-555 


11 


0, 


,57 


0. 


,72 


7.8 


2.2 


5.8 


6.7 


1.5 


631C 


B 


35.8±1, 


.6 


150-450 


5 


0, 


,30 


0. 


,68 


4.5 


5.3 


(16.1) 


11.9 


3.3 



34.8±0.9 4.69 



40.0±4.3 5.69 



37.2±2.7 5.18 



55.0±6.9 9.30 



51.8±6.8 8.62 



30.5±3.8 4.95 



36. 



39 



6±1.6 5.89 
7±5.7 7.18 



25, 



20, 



61 



Oil. 7 4.23 
1±1.7 2.78 
0±3.6 9.47 



37.7±2.4 5.76 



Grade: classification based on the number of a posteriori criteria checked (A all, B except at least one evidenced with paren- 
thesis); Fe±aFe: individual palaeofield strength determination (in fiT) and standard error associated jPrevot et al., 1985} ; 
AT: temperature interval of determination; n: number of consecutive points used for the determination; /, g and q: fraction of 
NRM, gap factor and quality factor, respectively jCoe et al., 1978l|Prevot et al., 1985} ; MAD: maximum angular deviation; 
a: angle between ChRM and origin; f}: angle between ChRMs of the sample and its sister specimen jPlenier et al., 2002} ; 



